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The influence of substitutional disorder on the transport properties of heavy-fermion systems is investigated. 
We extend the dynamical mean-field theory treatment of the periodic Anderson model (PAM) to a coherent- 
potential approximation for disordered strongly correlated electron systems. Considering two distinct local 
environments of a binary alloy AfSi-e with arbitrary concentration c, we explore two types of disorder: on 
the / site and on the ligand sites. We calculate the spectral functions and self-energies for the disordered 
PAM as well as the temperature dependence of the resistivity and the thermoelectric power. The characteristic 
concentration dependence as well as the order of magnitude of transport properties are reproduced for metallic 
heavy-fermion systems and Kondo insulators. In particular, sign changes of the Seebeck coefficient as function 
of temperature and concentration are observed. 
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PACS numbers: 71.10.Fd, 71.27.+a, 72.10.-d, 72.15.-v 



I. INTRODUCTION 



Heavy-fermion systems (HFSs) are characterized by quasi- 
particles with a very large effective mass at low temperatures 
T. This behavior occurs in many lanthanide (rare-earth) and 
actinide compounds and has its origin in the hybridization 
between the conduction bands and the local moments of in- 
completely filled / shells of the lanthanide or actinide ions.''^ 
These / electrons contribute substantially to the formation of 
the heavy quasiparticles. Measurements of the transport co- 
efficients on HFSs reveal characteristic anomaUes at low tem- 
peratures. The temperature dependence of the resistivity p (T) 
of metallic HFSs (such as CePds, CeAl3, CeCug, CeCu2Si2, 
UPd2Al3, etc.; see Refs. 2,3,4,5,6,7) shows a rapid increase 
with increasing T starting from a small residual value, which 
can often be fitted by a T^ law. The resistivity saturates 
at Tmax and decreases with increasing T for T > Tmax as 
p{T) ^ iog{Tmax/T), characteristic of the Kondo effect. In 
Kondo insulators,^ •^''^ however, a narrow gap opens''' at the 
Fermi energy at low temperatures. In such materials, e.g., 
SmBg or Ce3Bi4Pt3, the resistivity shows an activation be- 
havior for T ^0. Another interesting transport quantity is the 
thermoelectric power or Seebeck coefficient ^(r). At tem- 
peratures comparable to Tmax, the thermoelectric power can 
reach absolut values of about 50 /iV/K, often accompanied 
with sign changes at intermediate temperatures.''^ In Kondo 
insulators, even larger values of the thermoelectric power up 
to 300 /xV/K have been observed,"''^ which may be interest- 
ing for low-temperature thermoelectric cooling. 

In this paper, we examine alloys of heavy-fermion materi- 
als with arbitrary concentrations. Disorder is introduced by 
substitution of ligand ions {disorder on the ligand sites) or of 
lanthanide or actinide ions {disorder on the f site). A special 
case of disorder on the / site is given by replacing /-electron 
ions (e.g., Ce) by a certain amount of nonmagnetic impurities 
(e.g.. La) referred to as "Kondo holes." While we focus on 
local disorder only, band disorder in the Kondo lattice model 
has been recently considered'^ using a slave-boson mean-field 



approach. 

Disorder has strong impact on the transport properties 
such as p{T) and 5(7"). In substitutional alloys such 
as Lai_^Ce;,.Pd3,3 Ce.,Lai_.,Cu6,'' Ce.,Lai_.,Cu2.o5Si2,'''''^ 
Ui_vThvPd2Al3,'^, etc., the residual resistivity p(0) rapidly 
increases with increasing concentration x of the nonmagnetic 
impurities, and one obtains a crossover from the metallic 
p (r) behavior with a maximum to a monotonic curve, where 
p{T) decreases with increasing T. The thermoelectric power 
^(r) for CeCu2.o5Si2 shows a crossover from one case with 
a negative minimum, a sign change, and a positive maxi- 
mum to a behavior with two positive maxima and no sign 
change when substituting Ce by La.'"*-'^ For materials such as 
Ce3CuTPt3_vSb4, ligand alloying introduces a transition from 
a Kondo insulator to a dirty metal. '^ 

The basic model for a description of the electronic proper- 
ties of HFSs is the periodic Anderson model'^ (PAM). In a re- 
cent paper, '^ we studied the PAM within the dynamical mean- 
field theory '^'^^ (DMFT) and showed that the characteristic 
transport properties of HFSs can be understood within this 
framework. Within the DMFT, the lattice model is mapped on 
an effective single-impurity Anderson model^' (SIAM) by a 
self-consistency condition. We used the numerical renormal- 
ization group^^-^^ (NRG) as impurity solver for this effective 
SIAM. The NRG is a nonperturbative method applicable at 
very low temperatures which reproduces the correct charac- 
teristic low-temperature scale (Kondo temperature).^^ 

In the present paper, we extend the DMFT-NRG treat- 
ment to the disordered PAM. The disorder is studied within 
the coherent-potential approximation^^-^^ (CPA). The CPA 
was originally developed for alloys modeled by noninteract- 
ing particles. In the CPA, lattice coherence is restored by 
introducing an average potential. Janis and Vollhardt have 
pointed out that the local nature of the potential allows one 
to embed the CPA into the DMFT framework^^ emphasizing 
that the CPA is the best possible single-site approximation^^ 
for disorder. On the other hand, the DMFT can be inter- 
preted as CPA for periodic correlated electron systems as 
the DMFT reduces to the CPA equations in the absence of 



electron-electron interactions.^^ This allows for generaliza- 
tions of the CPA to disordered systems with finite Coulomb 
repulsion which has been exploited in the context of the Hub- 
bard modelP-^'^-^^'^^'^^ 

There exists a number of previous works in which the 

influence of disorder on the electronic properties of HFSs 
has been studied.33-34,35,36,37,38,39,40,4i,42 j^^^j ^f ^^^^^ p^. 

pers consider the case of disorder by Kondo holes only. 
The earliest work^^'^'* applied the alloy analog approximation 
to the PAM as a CPA application. Other treatments^^-^^-^** 
used the slave-boson mean-field approximation to the PAM, 
which enabled a straightforward application of the standard 
CPA.^^'^^ Schlottmann^'^ studied the development of impu- 
rity bands within the hybridization gap of Kondo insulators in 
the low-concentration limit using the iterative pertubation the- 
ory (ITP). Mutou'*'''*^ and Wermbter et aZ.^" calculated trans- 
port properties within the framework of the ITP and second- 
order perturbation theory, respectively. Recently, the low- 
temperature behavior of paramagnetic Kondo lattices upon 
random depletion of the local / moments has been investi- 
gated by Kaul and Vojta.^^ 

We extend the DMFT in order to describe a disordered 
PAM, based on the CPA for noninteracting systems in the 
spirit of Janis and Vollhardt.^^ In contrast to the standard CPA 
with static potentials, we use dynamical local potentials which 
include correlation effects. For reasons of simplicity, we call 
this extension "CPA", too, and stick to the name "DMFT" for 
pure systems without disorder. Here, we apply such a gener- 
alized CPA to the disordered PAM and calculate the transport 
quantities p{T) and 5(7") for disorder on the / site as well 
as for disorder on the ligand sites. We treat the whole range 
of impurity concentrations c E [0; 1] for different choices of 
the PAM parameters, starting with either a metallic HFS or a 
Kondo insulator for c = 0. We find that this CPA-NRG treat- 
ment of the PAM is able to explain the strong disorder depen- 
dence of the transport quantities of HFSs. 

A generalization of the CPA for binary alloys to infinitely 
many local environments has been considered in the context 
of ligand disorder of UCu5_rPd, by Dobrosavljevic and co- 
workers. '^^'^^'"'^^ Those locally different environments yield a 
distribution of hybridization strengths between the actinide 
and the conduction band at each site and, as a consequence, a 
distribution of Kondo scales. By averaging over these differ- 
ent low-energy scales, a non-Fermi-liquid behavior of trans- 
port properties was obtained. ^^'^^'^^ 

HFSs often exhibit symmetry broken phases such as 
superconductivity'*^ or antiferromagnetism,^ which are cer- 
tainly not yet properly understood. Such phase transitions 
at finite temperatures can be qualitatively understood within 
the PAM,'*'''50.5i.52,53,54,55,56,57,58 ^^^ jj^g complete phase dia- 
gram of the model in three dimensions has not yet been de- 
termined. A vanishing transition temperature as a function 
of an external control parameter such as magnetic field, pres- 
sure, or doping defines the quantum critical point (QCP), at 
which the transition is governed by quantum rather than ther- 
mal fluctuations. This sparked a lot of experimental ' as well 
as theoretical^^'^''"^' interest since it is believed that these fluc- 
tuations strongly affect the transport properties at a finite tem- 



perature in a cone close to the QCP^"'^' A proper microscopic 
description of such QCP requires the treatment of one- and 
two-particle properties on equal footing, which is beyond the 
scope of the present investigation and must be left for future 
research. Therefore, we restrict ourselves to the influence of 
disorder on transport properties in the paramagnetic phase of 
the model. 

The paper is organized as follows; In Sec. II, we introduce 
the model and the notations as well as our calculation meth- 
ods, the DMFT, the CPA, and the NRG as impurity solver 
Section III is devoted to the special case of Kondo holes: in 
Sec. Ill A, simplifications of the CPA equations are consid- 
ered, and in Sec. IIIB, we discuss the temperature-dependent 
spectral functions and self-energy of the PAM. These single- 
particle properties determine the transport properties via the 
Kubo formulas introduced in Sec. IV. We explicitly state 
equations for the resistivity p and the thermoelectric power 
S. In Sec. V, results on the temperature dependence of p{T) 
and 5(7") are reported, both for disorder on the / site and dis- 
order on the ligand sites. We conclude with a summary and 
outlook in Sec. VI. 



n. THEORY 



The essential features of heavy-fermion compounds are 
based on the interplay of localized, strongly correlated / elec- 
trons with broad conduction bands. At high temperature, 
the weakly coupled / electrons cause mainly incoherent and, 
with decreasing temperature, logarithmically growing spin- 
flip scattering for the conduction electrons. Below a char- 
acteristic temperature scale, a crossover to a coherent low- 
temperature phase is observed. The / electrons contribute sig- 
nificantly to the formation of heavy quasiparticles, while their 
moments are dynamically screened. The PAM takes these in- 
gredients into account, comprising of spin degenerate conduc- 
tion electrons, a lattice of correlated localized / electrons, and 
a hybridization (cf. Sec. II A). 

In this work, we investigate the influence of disorder in 
HFSs. In general, alloying with different ingredients A, B, 
C, . . . , with concentrations ca, cb, cc, . . . , respectively, de- 
stroys the lattice periodicity. Here, we consider only substi- 
tutional alloys of the type AcBi-c with concentrations ca = c 
and CB ~ I —c. 



A. Model 



This section is divided into two parts: In the first part, we 
will introduce the model for periodic systems, which is ex- 
tended to binary alloys in the second part. 



1. No disorder. 

The Hamiltonian of the simplest version of the PAM is 

given by'^ 



are possible with a probability c that is given by the relative 
concentration of subsystem A: 



^ = E ^la^U^la + Y L«fa«f-CT 



ka 



+1(4. 4) ^(t) with ^^(^^ 



(1) 



Here, c^^ (cl ) destroys (creates) a conduction electron with 

spin (7, momentum k, and energy Cj^ + Eca, where £ca in- 
dicates the band center. The energy e^^ denotes the spin- 
dependent single-particle /-level energy at lattice site /, h\^ = 
fjfjfia is the /-electron occupation operator (per site and 
spin), fia ififj) destroys (creates) an / electron with spin a at 
site i, and U denotes the on-site Coulomb repulsion between 
two / electrons on the same site /. The uncorrected conduc- 
tion electrons hybridize locally with the / electrons via the 
matrix element V^. 

Even though only a single effective / level is consid- 
ered, this model is quite general. It describes any heavy- 
fermion system with odd ground-state filling of the / shell, 
for which in a strong crystal field environment the degenerate 
Hund's rule ground state may be reduced to an effective spin- 
degenerate Kramers doublet. In addition, charge fluctuations 
to even / fillings leave the / shell in crystal field singlets. The 
Hamiltonian contains four energy scales. The interplay be- 



tween Eff, and U controls the average / filling as well as the 
local moment formation for large U and negative e^^. The 
Anderson width Fq = y^;rp()(0) determines the charge fluctu- 
ation scale of the / electrons, with po(0) being the density of 
states of the noninteracting conduction band of width D at its 
band center 

The total filling per site, «tot = L<j((«;a> + («/(7>)' is 
kept constant by a temperature-dependent chemical potential 
/i(r). We absorb the energy shifts into the band center Eca of 
the conduction band, as well as the / level Eff,- For «tot = 2 
and f/ = 0, the uncorrected system is an insulator at T = 0, 
since the lower of the two hybridized bands is completely 
filled. According to Luttinger's theorem, a finite U of arbitrary 
strength does not change the Fermi volume, which includes 
the full first Brillouin zone. As long as the ground state does 
not change the symmetry due to a phase transition, the system 
remains an insulator at arbitrarily large Coulomb repulsion. 
Therefore, the nonmetallic ground state of Kondo insulators is 
not correlation induced, but it is already present for the nonin- 
teracting system and is a consequence of Luttinger's theorem. 
For nonintegral values of «tot, the paramagnetic phase of the 
system must be metallic. 



2. Binary alloy. 

The simplest version (1) of the PAM can be extended to 
describe two subsystems by replacing the matrix '^ by a ran- 
dom potential %„. At each lattice site /, two different values 



^o 






with probability c 
with probability 1 - 



(2) 



In the case of no doping (i.e., 'f^ — Y^ or c = or c = 1), 
the extended model reduces to Eq. (1). 

While Eq. (1) includes possible Zeeman splitting of the en- 
ergies in an external magnetic field H, we set // = through- 
out the remainder of the paper and treat all properties as spin 
degenerate. In particular, we drop the spin index a from now 
on for the sake of simplicity. 



B. Dynamical mean-field tlieory 

Setting aside exact solutions in one dimension^^ using the 
Bethe ansatz for the Kondo lattice model, to our knowledge 
no exact analytical solution has been found for the model (1) 
with finite U . Therefore, one has to rely on suitable approx- 
imations for the PAM. An obvious first approximation is the 
assumption of a purely local, site-diagonal (i.e., k indepen- 
dent) self-energy, which for the PAM is even better justified 
than for other lattice models of correlated electron systems, 
as the first corrections are at least of order V^. Within a lo- 
cal self-energy approximation, the complicated lattice model 
can be mapped on an effective SIAM.^' Such a mapping was 
first used about 20 years ago in connection with the appli- 
cations of the noncrossing approximation^^'^-''^^'^^ (NCA) to 
the PAM.'''''^^'^^-'''^'^'' The effective site is viewed as a cor- 
related atomic problem within a time-dependent (or energy- 
dependent) external fleld^^'^^'^"'^' which is determined self- 
consistently, thus accounting for the feedback of the electron 
propagation through the lattice. ^^'^^'^^ Using the scaling of the 
tight-binding hopping parameter in large dimensions by Met- 
zner and Vollhardt,^^ Miiller-Hartmann^^ has proven that the 
local self-energy approximation becomes exact in the limit of 
infinite spatial dimensions. Therefore, the self-consistency 
condition of the DMFT,'^'^"'^^'^^'^^ which neglects spatial 
fluctuations in the single-particle self-energy, becomes exact 
in the limit d ^ °o. Within weak-coupling [/-perturbation the- 
ory, it could be shown^^'^^ that a local, ^-independent self- 
energy is a good approximation for realistic dimension d = 3 
as corrections due to intersite contributions to the self-energy 
are negligibly small. As a consequence, phase transitions re- 
main mean-field-like in DMFT since fc-dependent fluctuations 
are not included in a local approximation. ^"'^^'^^ 

The following exact relations for the conduction electron 
Green function Gc{k,z) and the /-electron Green function 
Gf(k,z) can be obtained for the PAM (1): 



Gik.z) 



Gff(k,z) Gfc(k,z) 
Gcf{k,z) Gcc{k,z)^ 







e/ + E/(^,z) V 

V Ec 



-1 



(3) 



where z is any complex energy off the real axis. Within a 
local approximation such as the DMFT, the fc-dependent /- 
electron self-energy 'Lf{k,z) is replaced by a ^-independent 
r/(z). From Eq. (3), one defines a self-energy of the conduc- 
tion electrons via 



Mz) = 



V' 



z-Sf-^fiz) 



which can include a simple k dependence through the hy- 
bridization matrix elements V^, here taken as a constant. 
For such a local self-energy E^ , the site-diagonal conduction- 
electron Green function Gcc can be written as a Hilbert trans- 
formation 

Gcciz) = ^£G,,(fc,z) =D[z-e,-Z,(z)], (5) 

k 

defined for arbitrary complex argument z as 

,Po(e) 



D{z) 



de- 



(6) 



where Po{(o) is the density of states of the noninteracting con- 
duction electrons. 

The DMFT self-consistency condition states that the site- 
diagonal matrix element of the /-electron Green function of 

~<\ocf 



the PAM must be equal to G'^J (z) of an effective site problem 



Gffiz)-^LGffiU)=G'rfi^), 

k 

1 



Gfn^) = 



(7) 



(8) 



with the same local /-electron self-energy T.f{z) for the lat- 
tice and the effective site. This defines the self-consistency 
condition for the functions E/(z) and A(z). 

We can put the site-diagonal Green functions into a Green 
function matrix 



G{z) = 



Gffiz) Gfciz) 
G,f{z) Gcciz) 



(9a) 



but it is sufficient to know one of the components because they 
are connected to each other: 



Gfc{z)=Gcf{z) = ^^Gcc{z) 



Hz) 



V 



V 



Gffiz). (9b) 



Furthermore, the self-consistency condition (7) extends to the 
full matrix, G{z) — G^°^{z)- It is possible to formulate the 
DMFT equation for the conduction-electron Green function 
instead of the /-electron Green function: 



Hz) 



-ef-T.f{z)-Gff{z)- 



V- 



j:.(z) + g,,(z)-i- 

This condition can also be written in the following form: 



(10) 



z 
z 



e/ + E/(z) V 
V e. 



G-\z) 





z~ec- 






Given the Green functions Gff{z) and G« (z), their spectral 
functions 



Pf{(o) = lmGff{co - /0+)/;r, 
Pc{o)) = lmGcc{oi - /0+)/;r 



(^^ determine the local occupation numbers 



nfic 



Y, d(of{oi-ii)pfic{(o), 



(12a) 
(12b) 



(13) 



where the spin sum contributes only a factor of 2 in the ab- 
sence of a magnetic field and /(ft)) denotes the Fermi func- 
tion. Then, the total filling per site is given by ntot — nf + n^. 
Particle-hole symmetry is reached at «tot = 2 and £/ — e^ = 
—U/2 for a symmetric po(ci)). As a matter of convenience, 
we will perform an integral transformation such that /i is ab- 
sorbed into Ef and the band center e^ ; all energies will be mea- 
sured with respect to /i. For a given lattice filling ntot, we have 
to adjust pL in addition to fulfill Eq. (7). 

Before we discuss the solution of the effective site, let us 
briefly comment on the implications of the analytical form 
of the conduction-electron self-energy (4). For a Fermi liq- 
uid at r ^ 0, the imaginary part of E^ vanishes quadratically 
close to the chemical potential, i.e., ImE/^(a) — /0+) °^ co". 
At particle-hole symmetry, E^ (z) diverges as 1 /z, leading to 
an insulator. Away from particle-hole symmetry, the denom- 
inator remains finite and the imaginary part also must have 
Fermi-liquid properties ImEc(a) — /0+) °^ ImE^(ft) — /0+) °^ 
(0^. The real part is very large, and therefore, the spectral 
function pdco) as well as Pf{(o) sample the high-energy band 
edges of po(co), yielding a hybridization gap. These analytic 
properties must be fulfilled by any approximate solution of the 
DMFT self-consistency condition (7) consistent with Fermi- 
liquid theory. 



C. Coherent-potential approximation 

The DMFT is only applicable to periodic systems. When 
we consider alloys or doped systems where the lattice pe- 
riodicity is destroyed, then the DMFT has to be extended. 
There are several possible approximations^^ for disordered 
alloys, e.g., the virtual-crystal method, the average t-matrix 
approximation,^** or the CPA.^^-^''-^^'^^ The CPA is the best 
single-site approximation^^'^** and can be seen as an extension 
of the DMFT for two subsystems. If there is effectively only 
one subsystem (as in the case of 1^ = % or c = or c = 1), 
the CPA reduces to the DMFT. 

The standard CPA self-consistency condition for noncorre- 
lated systems is given by 



ch-m-^'^^^w] m-^ 



CPAi 
1 



(11) 



(l-c){l-[rB-E^P^]^J [%-E^P^]=0, (14) 



where TF^^{z) is the self-energy (matrix) of the configura- 
tionally averaged Green function (matrix) 

%(Z) ^fciz) 
%f{z) %c{z) 



k 



Z 



-ECp^(z) 



-1 



(15) 



[as in Eq. (5), the CPA Green function matrix can also be writ- 
ten as a Hilbert transformation by replacing the k sum with the 
integral D] . 

To include correlation effects in Eq. (14), the potential ma- 
trix 1^ /g of subsystems A and B, respectively, has to be re- 
placed by the sum 



'^A/B^ 



"A/B- 



Z^/\z) 




(16) 



where Y,/ (z) is the local /-electron self-energy of subsys- 
tems A and B, respectively. Equivalent to Eq. (14) with corre- 
lation is the following CPA equation: 



^(z)=cG^(z) + (l-c)G^(z), 



(17) 



where C^ and G^ are the local Green function matrices (9) of 
effective site problems for the two subsystems A and B, given 
by the effective media A^(z) and Ab(z), respectively: 



Cf/r'i^ 



-.A/B 



Al. 



A B . , X ^A ti, si 



iz) 



^^^gT{z) 



M/B 



w. 



V/ 



A/B 



Va/b 



-.A/B 



iz). 



(18a) 
(18b) 



Note that the determinant of ^ [Eq. (15)] can be written in 
the following simple form: 



det^(z) 



, _ vCPA ■ 



(19) 



determine the occupation numbers 



nf/c ^ 



d(of{co~^)pfJ^{(o). 



(23) 



In analogy to Eq. (4), we define 



-CPA 



(z) 






(24) 



The CPA self-consistency cycle consists of two parts: (I) 
After solving the two impurity subsystems A and B for the 
given effective media, the CPA self-energy matrix is deter- 
mined by Eqs. (17) and (20). (II) A new estimate for the CPA 
Green function matrix is obtained by Eq. (15). Both matri- 
ces are required to calculate r(z) — and thus two new effective 

A/B 

media — by Eq. (20). The cycle is completed by shifting e/ 

and Ec such that N[G] = «tot and jU = 0. 

Although the CPA includes the DMFT as end points for 
the concentrated systems (c = and c ~ I), conceptual dif- 
ferences in disordered systems arise from the occurrence of 
the additional quantities. In the CPA, we must distinguish be- 
tween the local matrix 7^ /g, which describes the properties of 
a site of type A/B and differs at different sites, and the con- 
figuration averaged A:-independent lattice self-energy lP^'^{z), 
which is equal for all lattice sites. Moreover, the configura- 
tionally averaged CPA Green function restores translational 
invariance in the lattice and is not equal to either local Green 
function G'^'^{z)- Even though the CPA yields only one con- 
figurationally averaged dynamical field r(z), the local dynam- 
ics of each site is determined by different media A^/3(z). Only 
for c = and c = 1, the CPA equation (17) is trivially fulfilled 
and TP^^{z) coincides with the local self-energy as well as 
the CPA Green function ^(z) coincides with the one of the 
effective site. 



Using this expression, it can be shown by direct calculation 
that 



z 
z 



-E^P^(z)-^-'(z) 





r(z) 



(20) 



has got one component only. Therefore, we obtain the same 
structure as in the case of the DMFT, Eq. (11). This leads 
to a DMFT-like self-consistency for the effective medium by 
setting 



^a/b{z) ■■-- 



V, 



A/B 



A/B r^, N ■ 

z-Ec' -r(z) 



(21) 



This ensures that for pure systems (c = or c = 1) the DMFT 
limit as given by Eq. (1 1) is recovered. 

As in Eqs. (12) and (13), the CPA spectral functions 



py'^ico) = Im%(co - iO+)/7t, 
p,cPA(co)=Im^,,(co-/0+)/7r 



(22a) 
(22b) 



D. Impurity solver 



A/B 

The local /-electron self-energy T.J entering the CPA 

equation (17) via the local Green function matrix G^'^ 
[Eq. (18)] is obtained for an effective site defined by the lo- 
cal Hamiltonian'^''^'^"-^^-^^'^^ H^^ for subsystems A and B, 
respectively. 



+ Y. de Vaib\Jp%\>^) {dUa + fUoe) ■ (25) 

The coupling of the / electron to a fictitious bath of "con- 
duction electrons" created by djjg is described by an energy- 
dependent hybridization function 

^yI/Bpt«'i£) - ImA^/g(e - /0+) . (26a) 



Using Eq. (21), we can write the density of states of the ficti- 
tious bath as 



PeT(e) = - Im[e - /0+ - e^'"" - r(e - /O^ 



(26b) 



Yl ^ 2t*YQ/V2Ti 



IQTl 



SMTl 



K ss 5.641 Q. In order to make 
contact to the experimentally relevant parameter regime, we 
assume To = 100 meV, which translates into a temperature 
scaleof ^wll60K. 



We accurately solve the Hamiltonian (25) using Wilson's 
NRG approach. ^^'^^ The key ingredient in the NRG is a log- 
arithmic discretization of the continuous bath, controlled by 
the parameter^^ A > 1. The Hamiltonian is mapped onto a 
semi-infinite chain, where the Mh link represents an exponen- 
tially decreasing energy scale D^r ^ A^'^'^. Using this hier- 
archy of scales, the sequence of finite-size Hamiltonians J^/f^ 
for the A^-site chain is solved iteratively, discarding the high- 
energy states at each step to maintain a manageable number of 
states. The reduced basis set of .^^ thus obtained is expected 
to faithfully describe the spectrum of the full Hamiltonian on 
the scale of Dn, corresponding^^ to a temperature 7)v ^ Dn 
from which all thermodynamic expectation values are calcu- 
lated. The energy-dependent hybridization function A(z) de- 
termines the coefficients of the semi-infinite chain. '^'^^ For 
further details, we refer to the recent NRG review by Bulla et 

The finite-temperature NRG spectral functions are calcu- 
lated using the recently developed algorithm from Refs. 82 
and 83. The usage of a complete basis set of the Wilson 
chain originally derived for real-time dynamics of quantum 
impurites out of equilibrium^"*-^^ ensures that the spectral sum 
rule is exactly fulfilled and that the NRG occupancy is accu- 
rately reproduced for arbitrary values of A and number of re- 
tained states Ns- The NRG spectrum is broadened in the usual 



way 



81,86,87 



by a Gaussian at higher frequency and a Lorentzian 



for low-frequency excitations \E \ 

. , i^e-IK'sC"/^)/^]' for|£|>L,T 
\li (.-EYHLrTY ior\E\<U.T. 

As long as not stated otherwise, we use 

b = 0.6, Lr = 1, K^ y^LTexp{-b^/4-)/b. (27b) 

A/B 

The local /-electron self-energy E^' is determined by the 
exact ratio 



^NRG^,^ 



(28) 



derived via equation of motion technique,**^ where the 
con-elation functions M^^°{z) = {{foflaf-c\fl)){z) and 
F^^'^{z) = {{fo\fl)){z) have been obtained from NRG spec- 
tral functions. 

We use a Gaussian model density of states 



Po(e) = 



exp[-(£A*)V2] 
t*V2n 



(29) 



for the unperturbed conduction-electron system. In the fol- 
lowing, we set \/2f* — IOFq and measure energies in units of 
To = ;rV'o po(0) which defines a reference hybridization Vo by 



III. KONDO HOLES 
A. Analytical considerations 

A special case for disorder on the / site is given by in- 
troducing Kondo holes, i.e., lattice sites without / electrons. 
We describe the two subsystems with identical parameters 
{V ■.= Va = Vb, Ec := e^ = ef , and therefore A := A^ = A^) 
except that the /-level energy of system B is shifted toward 
infinity. Additionally, the total electron density is reduced by 
one electron per B site to «tot = ntot(c=l) — (1 — c). 

Because there are no / electrons in system B, the 
conduction-electron self-energy Ef [Eq. (4)] vanishes. By 
Eqs. (9) and (11), we conclude that the Green function ma- 
trix of system B contains only the cc element: 



G'iz) 











A{z)/V' 



(30) 



As a consequence, the configurationally averaged /-electron 
Green function of Eq. (17) is given by ^//(z) = cG4^(z), van- 
ishing for c ^ 0, and the configurationally averaged band- 
electron Green function '^cc{z) interpolates between G^^.(z) 
and the free medium with band center at Ec- 

By inserting Eq. (17) into Eq. (20), we conclude that in the 
case of Kondo holes the CPA self-energy matrix reduces to a 
scalar self-energy 



-CPA 



(z) 



T.7Mz) V 



'ff 



ec 



(31a) 



with 



^f^iz) 



'(z)]// = z-Aa(z)- 



cGUz) 



(l-i)[z-A(z)] + i[4+E)(z)]. 



(31b) 



There are two contributions to the imaginary part of ES:?'^: 
a term ilmE4 and a term i-^ImA. If E4 remains Fermi- 

c I c J 

liquid-like for O), T ^ 0, a finite life time is introduced by the 
additional ^-^ ImA(O) for c < 1. We thus interpolate between 
the following DMFT limits: 



lim 



hm 



lim 
lim 



c^)nz) 



civru^) 



-// 



pCPA 



(z) 



^'\z) 



^ef + mz), 


(32a) 
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(32b) 


Gjfizy 


^^^iz), 


(32c) 


= V'-Gjf{z). 


(32d) 



Because the CPA mixes two subsystems, it is not clear 
how to define a low-temperature scale for the complete sys- 
tem. Each effective site has its own low-temperature scale 
that characterizes the local dynamics. In the case of Kondo 
holes, we are in a better situation. As stated above, the con- 
figurationally averaged /-electron Green function is directly 
connected to the /-electron Green function of subsystem A: 
^/■/(z) = cG4,(z), there are no / electrons in system B and, 
hence, no special low-temperature effects. Then we can take 
the low-temperature scale of subsystem A to define a low- 
temperature scale Tiow for the complete system which is only 
defined up to a constant factor We will use the renormal- 
ization of the Anderson width by the quasiparticle spectral 
weight 



7b = ro 



1- 



dReI.j{(0) 



do) 



m,T^O 



(33) 



as our choice of such a low-temperature scale 7b o= Tiow for our 
numerical analysis. ^^ It is related to the mass enhancement 

m*/m = ro/To. 



B. Spectral properties 

In this section, we take a look at the spectral properties of a 
system with Kondo holes. 

The typical structure of the /-electron spectral function 
^p9^^(a)) = pf{co) is shown in Fig. 1. A pronounced peak 
structure dominates the low-energy part of the spectrum in the 
vicinity of the chemical potential. In addition, we observe 
two shallow high-energy peaks: one at Ef below ji, and one 
at Ef + U, which corresponds to double occupancy of the / 
levels. It is always correctly positioned by the NRG, indepen- 
dent of the value of U, but with a linewidth too large due to 
the NRG broadening procedure of Eq. (27) (see Refs. 87 and 
88 for details). 

The dependence on the concentration c is small, but the 
peak width is enlarged with decreasing c. This indicates an 
increasing low-temperature scale; indeed, the scale given by 
Eq. (33) and shown in the inset of Fig. 1(b) increases with 
decreasing c. 

The band-electron spectral function is depicted in Fig. 2. 
Equivalent to the vanishing of / electrons from the system 
with decreasing c is the disappearing of the gap structure and 
the crossover to a free medium of Gaussian shape as discussed 
in the previous section. 

The behavior of the spectral functions is reflected in the 
conduction-electron self-energy -E^^'^(z) (cf. Fig. 3). At low 
temperatures and c ~ 1, we observe a characteristic (O^ de- 
pendence at the chemical potential jJ. ~Q, indicating a Fermi 
liquid. At c = 0, the self-energy is proportional to G'lx. The 
imaginary part of the conduction-electron self-energy at a 
temperature of 274 K can be seen in Fig. 3(b). Additionally, 
the derivative of the Fermi function, /', is shown. The prod- 
uct of/' and ImE^^'^ enters the transport integrals in Sec. IV. 
That is, only the part of ImE^'''^ around O) = is relevant for 




FIG. 1: Spectral density Pf(co) for different concentrations of sys- 
tem A, c e {lO""*, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,1}, calcu- 
lated with CPA-NRG for Va = Vb = Vq, U /Tq = 10, ej-£^ = 

—U/2, e? = oo, chemical potential /x = 0, and a filling ntot = 
1.6 — (1 — c) at a small temperature 0.2 K. (a) Overview of the 
generic three-peak structure, (b) Peak at O) = in more detail with 
corresponding temperature scale To [Eq. (33)] in the inset. NRG pa- 
rameters: number of retained NRG states, A^, = 800, A = 1.6, and 

8/ro = io-\ 
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FIG. 2: Spectraldensitypc(co), with all parameters as in Fig. 1. With 
c -^ 0, the gap at tt) = vanishes and we reach a Gaussian curve. 
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FIG. 3: Band self-energy \'^^^ for different concentrations of sys- 
tem A, ce {10-4,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1}, and all 
parameters as in Fig. 1; (a) at T ~ 0.2 K, (b) at T ~ 274 K, with 
additional derivative of Fermi function, /', showing the relevant part 
of the energy axis. 



transport calculations, and in this range, the imaginary part of 
the conduction-electron self-energy is nearly independent of 
the concentration c. 



IV. TRANSPORT THEORY 

To describe the electronic transport within the PAM, we 
start from the standard relations^" for the generalized transport 
coefficients, according to which the electrical current density 
J and the heat current density q depend linearly on the electric 
field E and the temperature gradient VT: 



where Vq = cr" is the volume of the unit cell and A^ counts the 
number of lattice sites. It has been shown-''* that the current 
operator of the PAM has two contributions: a conduction- 
electron part and a part proportional to Wr. The (bare) / 
electrons do not appear in the current, since they do not dis- 
perse. For a ^-independent hybridization, only the conduction 
electrons carry the electrical and heat currents as follows: 



to 






(36) 



^V^e^ is the group velocity. Hence, the current- 



where vr 

K II K K 

susceptibility tensor ((j\j')){z) 
hole Green function 



is connected to the particle- 



(01/)) w 



aa'kk' 



vvv- ((cI CJ \c 



I'L'k'a'Mz). (37) 



In a cubic crystal, the conductivity is isotropic: (J„ft(a)) = 
(t((b)]1. From now on, we will consider only the xx compo- 
nent of the conductivity g{(o) = Gxxi(o)- 

In general, the full two-particle Green function 
((cl cr let c.-, ,))(z) involves vertex corrections which 

^^ to kal 1,1 (ji k'a"i'-~' 

reflect residual particle-particle interactions.^*' However, in 
the limit d ^ oa, \i was shown that current operator vertex 
corrections vanish. ^"^'^^ Thus, it is consistent with the DMFT 
assumption of a ^-independent self-energy that these vertex 
corrections vanish for any lattice model of correlated electron 
systems. For the special case of a local approximation for 
the PAM, this was already shown in Refs. 96, 97, and 78, as 
for symmetry reasons ^jvj|VjpGj(z + a))Gj(z) = 0. The 
CPA maintains the locality of the self-energy and we can 
neglect vertex corrections in the case of the CPA, too (see 
also Ref. 79). All DMFT formulas can be used for the CPA; 
we only have to substitute the DMFT self-energies with CPA 
self-energies. In both cases, we obtain 



((au1))(«+/o-)=^i:© 



dk. 



d(o'f{co')p,{er,(o') 



X [Gcc(k,a' + oi + iQ+) 
+ G„(^,co'-fi)-/0+)]. 



(38) 



Within the DMFT or the CPA, the lattice one-particle Green 
function depends only on the (complex) energy z and bare 
band dispersion e^: Gcc{k,z) — Gcc{£i,z)- Then 



J = LuE 



-Lu{-^yT) 



q = L2iE+L22{-Y'VT) 



(34a) 
(34b) 



All coefficients are calculated within the linear response ap- 
proach, starting from similar Kubo formulas. ^*''^' For symme- 
try reasons, Ln = L21 must hold. 

For example, the real part of the frequency-dependent (op- 
tical) conductivity tensor^**'^"'^^ c('^) = Lii{(o) is related to 
the current-current correlation function and is written as 



<yapio)) = - 



1 



CONVo 



Im((j„|yb)(« + /0+), (35) 



with 



^?(tl ^<^'' 



depoie)A{£), (39) 



A,w4s(5)"'"'-'!) 



(40) 



We use a Gaussian model density of states [Eq. (29)] for the 
unperturbed conduction-electron system which is appropri- 
ate for a (i-dimensional hypercubic lattice in the limit d -^ 



ooj-^ and po(e) has been evaluated approximately in large 
dimensions'^ as 






(41) 



on a hypercubic lattice. Then, Eq. (38) can be reduced to a 
sum of Hilbert transforms, which is defined in Eq. (6). 

By taking the limit © ^ 0, the static conductivity O — Ln 
and the thermoelectric power 



S 



1 Ln 



e kgTLi 



(42) 



are obtained.'"''^ The thermoelectric power S is defined as the 
proportionality constant between an applied temperature gra- 
dient and the measured voltage drop in the absence of a cur- 
rent flow. The Peltier coefficient, given by the ratio of heat 
and electrical current, is related to the thermoelectric power 
by n = TS. 

In the limit o ^ 0, we get for the generalized transport 
coefficients the following: 



(cT=)Ln = 
Ln = 

L22 = 



1 
ha 

1 
ha 

1 
ha 



[-/{(o)] e- T(a)) da, (43a) 

h/'(a))] e(S) t{co) da, (43b) 

[-/'(«)] 0)2 t(«) da. (43c) 



Here, /' is the derivative of the Fermi function; thus, we have 
the limits 

2 
Lii(r=0) = ^T(0), Li2(r=0)=0, L22(7'=0)=0, 

(43d) 



ha 



and T:{a) represents a generalized relaxation time defined as 

Po(e)p,"(e,«ye 
(f*)2 \dlmD{z) lmD{z) 



m-^ivf 



%d 
which is to be evaluated at ; 
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(44) 



(DMFT) 



fi) + /0+-e,.-E,.(fi) + /0+) 

a + /0+ - E^P'^(fi) + /0+) - EP^(« + /0+) (CPA) 



Here, e^ corresponds to the component E^^^'^, and Ec and E^^'^ 
are defined in Eqs. (4) and (24), respectively. For small y = 
Imz, e.g., in the Fermi-liquid regime, the approximation 



T(a)) 



{f) 



*\2 



Pq{x) p'^{x) 



L y 



y+o{y^) 



(45) 



is valid. In lowest order, we have T{a) °^ l/ImE^ (ct) + /0+) 
and the linearized Boltzmann transport theory is recovered. 
Though we are using the full form of Eq. (44) for calculations, 
it is sufficient to consider the lowest order only for discussion 
of the properties of the transport coefficients because it is the 
main contribution. 



V. RESULTS FOR THE TRANSPORT PROPERTIES 

All transport calculations rely on the results for the single- 
particle Green functions of the PAM. As stated in Sec. II D, 
we use a Gaussian model density of states [Eq. (29) with 
^/2t* = lOFo] for the unperturbed conduction-electron sys- 
tem and measure energies in units of Fq = 100 meV. In 
particular, this defines a reference hybridization Vq by Vq = 
2f *Fo/V27r = lOFg/ v^ « 5.64F2. 

If we assume one electron per unit cell of the volume a^ 
(a = 10^'" m), the resistivity p = a^' has the natural unit 



Oq ' = ha/e~ w 41 iiQ. cm. 



(46) 



Note that eL\2/kBTL\\ is dimensionless and kg/e k, 
—86 jUV/K. Therefore, the thermoelectric power is given in 
absolute units; only the scale of the temperature axis must be 
fixed by experiment. With our choice of Fq, it is given in units 
of ^w 1160 K. 



A. Disorder on the / site 

In this section, we restrict ourselves to the special case of 
Kondo holes, as introduced in Sec. III. In particular, the /- 
level energy of system B is shifted to infinity. The hybridiza- 
tion y = y^ = Vfi is set to the reference value Vq- 



1. Resistivity 

In Fig. 4, the resistivity is displayed for U /Yq = 10, Ex — 

e^ = — C//2, and various values of the concentration c of sys- 
tem A. The filling is set to «tot = 1.6— (1 —c). 

In the pure case, c = 1, we reproduce the typical behavior 
of metallic heavy-fermion systems^'^'^'^ within ourCPA-NRG 
treatment: a resistivity increasing with increasing T for low T, 
a maximum of the order of 100 ptD. cm at a characteristic tem- 
perature Tmax, and a p{T) (logarithmically) decreasing with 
increasing T for T > Tmax- Note that even for c = 1, a finite 
p (r = 0) is obtained for T ^Q. This is due to the fact that we 
are using a finite 5 = IO^^Fq contributing to the self-energy 
imaginary part. In addition, the broadening of Eq. (27) and 
the limited accuracy of the NRG in the regime a <T further 
enhance the self-energy imaginary part; for a discussion, see 
also Ref. 18. 

In the disordered case, c < 1, we distinguish the two 
regimes T < 7'max(c=l) and T > rmax(c=l) to obtain a quali- 
tative classification for all concentrations c. At high tempera- 
tures T > Traax, the logarithmic gradient of p(r) has its origin 
in the weak incoherent scattering of conduction electrons by 
the / moments. With increasing concentration 1 — c of Kondo 
holes, there are less scattering events. Indeed, the resistivity 
decreases and scales linearly with the concentration c as de- 
picted in Fig. 4(b). The c dependence of the resistivity p as 



well as of p/c at a fixed temperature T 
plotted in greater detail in Fig. 5(b). 
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FIG. 4: (Color online) Resistivity p{T) as a function of T for dif- 
ferent concentrations c of system A, calculated with CPA-NRG; all 
parameters as in Fig. 1, i.e., U /Tq = 10, ej — £^ = —U/2, e? = cc, 
and Mtot = 1.6— (1 — c). In (b), the same data are plotted as p{T)/c. 
The dash-dotted lines mark two cuts through the resistivity curves 
which are shown in Fig. 5 in more detail. 



This can be understood using the transport equation (43a) 
with abbreviation (46) as follows: 



- = C7, 



-1 



[-f'{C0)\cT{(0)d(0 



Because -\n\Lc{(o) is nearly independent of c for tempera- 
tures T > Tmax [cf. Fig. 3(b)], also ct{(o) does not depend on 
c [see Fig. 5(d)]. 

At temperatures T < Tmax, conduction electrons and / elec- 
trons tend to form heavy quasiparticles with Fermi-liquid 
properties, leading to a small resistivity for c = 1 and T ^0. 
The introduction of Kondo holes destroys the lattice periodic- 
ity and, therefore, the resistivity increases; for c < 0.5, the 
nonmonotonic course changes to a behavior which is typi- 
cal for a dilute distribution of impurity scattering centers [cf. 
Fig. 4(a)]. As a consequence, the resistivity at T < Tmax does 
not scale with the concentration c of Kondo holes. Instead, 
the normalized residual resistivity p{T=Q)/c depends on c 
and reaches a maximum at c = 0.5 [cf. Fig. 5(a)]; for small c 




FIG. 5: [(a) and (b)] Resistivity as a function of the concentration 
c at r ~ 0.2 K and T ~ 274 K, respectively: p{c) (dashed line), 
p{c)/c (full line), and approximation oc 1/(1 —0.6c) (dotted line), 
(c) Illustration of the proportionality p(T -^ 0)/c <x ImZ^^'^(0)/c 
[cf. Fig. 3]. (d) Generalized relaxation time fT(co) for T ~ 274 K. 
The additional Fermi derivation shows the relevant part of the energy 
axis. The temperatures correspond to dash-dotted lines in Fig. 4. 



it can be approximated by 



with rn := lim . 

l-0.6c c^Q c 



(47) 



This is consistent with other calculations^^'-'^-^^ and can be un- 
derstood using the CPA equation. The approximation Eq. (45) 
leads to p(r=0)/c = (7o"VcT(0) ex -ImE^P'^(0)/c [using the 
Umit (43d) with abbreviation (46)]. With the CPA equation 
(31), we obtain 



1 



lmi:^^^(0) 



1 



v^ 



c''"-lCPA(0) 



Im- 



V^GJfiO) 



1 



-cAa(0)G').(0)' 



(48) 

Because A^(0)G4(0) does not vanish, ImE^'''^(0)/c as well 
as p(T=0)/c acquire a 1/(1 — ac) dependence. 

A word is in order about the relation of the concentration 
dependence of the resistivity to the typical scaling behav- 
ior for heavy-fermion materials, e.g., in CevLai_vCu6.^ The 
reported magnetic resistivity contribution normalized onto 
the Ce concentration decreases monotonically with increas- 
ing cerium concentration.^ In Cej^Lai_;fCu2.o5Si2, however, 
a nonmonotonic concentration dependence of the magnetic 
part of the resistivity, normalized to its room temperature 
value, was found. *'*'^ Our calculations describe consistently 
the crossover from a concentrated lattice system to a di- 
lute distribution of impurity scattering centers. As discussed 
above, we can clearly pinpoint the origin of the maximum of 
the normalized residual resistivity p(0)/c at concentrations 
c = X = 0.5. It arises from multiple scattering processes at 
intermediate concentrations, and the CPA equation yields a 
low-temperature enhancement factor of about 1/(1 — ac) [see 
Eq. (47)]. In the dilute limit, p(0)/c becomes independent 



11 



500 



-^ c=l. 


r= 1 


-«-<: = 0.9, 


1- = 0.95 


-H c = 0.8. 


|. = 0.91 


-•- c = 0.7. 


I. = 0.86 


-♦- c = 0.6. 


i> = 0.82 


-•- c = 0.5. 


1. = 0.764 


^*- c = 0.4. 


1. = 0.73 


-•r- c = 0.3. 


1. = 0.69 


-*- c = 0.2, 


I. = 0.65 


-«- c = 0.1. 


I. = 0.61 


-»- c = lc-4 


,v = 0.574 




r[K] 



FIG. 6; (Color online) Resistivity as in Fig. 4(b), but with c- 
dependent hybridization V~ = v{c)Vq. 



of the concentration where the conduction band self-energy is 
proportional to the single-particle f-matrix.^^ 

A weak dependence of the low-energy scale as a function 
of doping has been reported in experiments.^'^ To simulate 
these findings, we also studied an explicitly concentration- 
dependent hybridization: V^ is scaled down in such a way that 
the same temperature scale Tq is obtained for all concentra- 
tions. This is achieved by a linear interpolation of V^ between 
the concentrated limit and the dilute limit. Since La increases 
the lattice spacing^'^^ of the materials, a decrease of the hy- 
bridization strength with increasing concentration 1 — c of the 
Kondo holes is consistent with the experimental findings. In 
Fig. 6, the resistivity is displayed for the same parameters as in 
Fig. 4, but for c-dependent hybridization. As a consequence, 
the calculated curves for the resistivity scale with cv{c). Since 
we do not know how to relate the explicit change of the hy- 
bridization to a change of the unit cell volume, we are not able 
to rescale the absolute unit. 

The experimental determination of the resistivity faces 
some severe limitations in determining the absolute value of 
magnetic contribution to the resistance. Setting aside any 
problems stemming from the subtraction of the phonon con- 
tributions by using a proper reference material,^''"* three major 
sources of errors remain: (i) the accuracy of the geometry fac- 
tor relating the measured resistance to the resistivity, (ii) the 
accuracy of the determination of the concentration c, and (iii) 
dealing with grain boundaries typical for such heavy-fermion 
alloys. In particular, the error of the concentration value c 
increases significantly in the dilute Umit, which has a pro- 
found impact on the absolute value of p{T)/c. In addition 
grain boundaries in samples might require the subtraction of 
the additional boundary resistivity in order to extract Pmag(7^) 
instead of normalizing the resistivity to its room temperature 
value, since grain boundaries are highly sample dependent. 
Therefore, we limit ourselves to the reproduction of the qual- 
itative features of the reported experimental behavior. 

We also investigated the influence of the local Coulomb re- 
pulsion on the transport properties. In Fig. 7, the resistivity is 
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FIG. 7: (Color online) Resistivity p (T) as a function of T for concen- 
tration c = 0.8 of system A, calculated with CPA-NRG for different 



(//Fo, e^ — e'^ = —U/2, e? = oo, chemical potential /x = 0, and a 



V 



displayed for a fixed concentration c = 



: 0.8, ej 



-U/2, 



filling 7itot = 1.4. NRG parameters: number of retained NRG states, 
Ns = 800, A = 1.6, 5/ro = lO"''', and L„, = (i.e., no Lorentzian 
broadening). 



and various values of U /To- The filling is set to «tot = 1-4 
to match the previous setting «tot = 1.6 — (1 — c) in the case 
of fixed U /Tq = 10. Here, we reduced the small imaginary 
part S/Fo to 10^"' to reduce the error in the residual resis- 
tivity to a minimum. While it turned out to be sufficient to 
use S/Tq = lO^-' for values of the Coulomb interaction up to 
U/To = 10, this is no more the case for U /Tq > 10. 

For a particle-hole symmetric conduction band, no even- 
odd oscillations in the NRG spectral function have been 
reported. ^^ On the contrary, we observe even-odd oscillations 
in the NRG spectral function of the effective site with respect 
to the number of NRG iterations, which is related to the tem- 
perature of the calculation in the usual way.^''^' This feature 
is related to the two strong coupling fixed points^^ of the NRG 
and can show up for particle-hole asymmetric band density of 
states. We used the standard technique for thermodynamic al 
properties such as entropy^^'^' and average the resistivity as 
a function of temperature for even and odd numbers of itera- 
tions. 

The temperature T^ax, at which the resistivity has its maxi- 
mum, is shifted to lower values for increasing U; one obtains 
an exponential dependency of Tmax on U. For sufficiently 
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FIG. 8: (Color online) Thermoelectric power S{T) as a function of T 
for fixed U /Tq = 10 and different concentrations c of system A; all 
parameters as in Fig. 4, among others «tot = 1.6— (1 — c). 



strong U, the peak height at Ty^ax is nearly independent of 
the value of U. This is characteristic for the disordered sta- 
ble moment regime, in which the / occupation reaches inte- 
ger valence, and the increase of U only reduces the effective 
Kondo coupling / tx V^/U. Charge fluctuations are strongly 
suppressed at low temperatures, and one reaches a universal 
regime for a fixed concentration. This is seen more clearly 
from Fig. 7(b), which shows the scaling properties plotting the 
resistivity p{T) versus T /Tq. We have defined Tq in Eq. (33) 
as an effective low-temperature scale of subsystem A in order 
to clarify the scaling properties of transport properties. We 
certainly do not imply universality of all low-energy proper- 
ties. We note, however, that Tq is of the order of the position 
Tmax of the maximum of the resistivity. While for U /Tq < 8, 
where high- and low-temperature scales are not very well sep- 
arated, the maximum of the resistivity and the peak height 
show a U dependence, and we reach a universality regime for 
large U. 



2. Thermoelectric power 

The thermoelectric power S{T) measures the ratio between 
electrical and heat current divided by the temperature, and 
its sign is related to the integrated particle-hole asymmetry 
relative to the chemical potential. In Fig. 8, S[T) is plot- 
ted for fixed U /Tq = 10, e4 — e^ = — C//2, and various val- 
ues of the concentration c of system A. The filling is kept 
at «tot = 1.6 — (1 — c) and the /-level energy of system B 
is shifted to infinity. This figure is accompanied by Fig. 9. 
Here, 5(7") is plotted for fixed filling ntot = 1-4, c = 0.8, 



e^ — e^ = —U/2, and various values of the Coulomb inter- 
action U [for the low-temperature scale Tq, see Eq. (33)]. 

We obtain very large absolute values for ^(r) (of the mag- 
nitude lOOjLiV/K; see also Ref. 18). Note that the thermoelec- 
tric power is obtained in absolute units, as already mentioned 
in the beginning of this section. Similar to the resistivity, the 
thermoelectric power exhibits a low-temperature peak which 
is correlated with the maximum of the resistivity, which is an 
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FIG. 9: (Color online) Thermoelectric power S{T) as a function of T 
for fixed concentration c = 0.8 of system A and different U /Tq; all 
parameters as in Fig. 7, among others ntot = 1.4. 



analytical consequence of Eq. (42). Therefore, the position 
of this low-temperature peak depends on the low-temperature 
scale Tq, which varies with the concentration c. This can be 
seen in Fig. 9(b), where the thermoelectric power is shown on 
a rescaled axis T /Tq. 

In addition, we observe a second extremum at a very high 
temperature independent of the concentration and the low- 
temperature scale which results from the charge fluctuations 
on the energy scale ej — jJ.. This maximum moves to higher 
temperatures with increasing U, as shown in Fig. 9. 



B. Disorder on the ligand sites 

In this section, we use the full matrix version of the CPA 
of Sec. lie, where the two subsystems may have different 
parameter matrices f^ and %; i.e., we introduce disorder on 
the ligand sites. We consider two cases: At first, we keep 
the same hybridization for both subsystems, Va =Vb = Vq; 
afterwards, we use a site-dependent hybridization Va 7^ Vb. In 
both cases, we fix the interaction atU ~ IOTq. 
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FIG. 10: (Color online) (a) Resistivity p{T) and (b) tliermoelectric 
power 5(r) as functions of T for different concentrations c of system 
A, calculated with CPA-NRG for Va = Vb = Vq, U /Tq = 10, e4 - 



•/ 



e;* = -U/2, e 
and a filling «tot 
NRG states, A'^ 



Fo, Er = £?, chemical potential ^ = 0, 



^R( 



= 2 — 0.4c. NRG parameters: number of retained 
800, A = 1.6, and S/Fq = 10"^ 



1. Resistivity and thermopower for Va = Vb = Vq 



FIG. 11: (Color online) (a) Resistivity p{T) and (b) thermoelectric 
power S{T) as functions of T for different concentrations c of sys- 
tem A, calculated with CPA-NRG for Va = Vb = Vq, U/Fq = 10, 
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tential 11 = 0, and a filling /itot = 1.8 — 0.2c. The inset of (b) shows 
the antisymmetrized relaxation time t(co) — t(— ft)) [Eq. (44)] in the 
vicinity of the chemical potential for the temperature Ti = 0.48 K. 
NRG parameters: number of retained NRG states, Ns = 800, A = 1.6, 
and5/ro = 10"^ 



We start with the investigation of the transport properties as 
a function of disorder in the case in which ligand substitution 
drives the system from metallic to Kondo-insulating behavior. 
The concentration-dependent filling is set to ntot = 2 — 0.4c, 
and the system evolves into a Kondo insulator by reducing 
c. In Fig. 10, the resistivity and the thermoelectric power are 
displayed for e4 — e^ = —U/2, e? — ef = — Fq, and various 
values of the concentration c of system A. The curves for c = 1 
are identical to the c = 1 curves in Figs. 4 and 8. We plot the 
resistivity on a log-log scale in Fig. 10(a) to cover the almost 
3 orders of magnitude of resistivity change. 

In contrast to the previous case of / disorder introduced by 
Kondo holes, the number of Kondo scatterers remains con- 
stant as a function of the concentration c. Therefore, the resis- 
tivity does not vanish for c ^ 0. For temperatures T > Tmax, 
the absolute values of the resistivity remain nearly the same; 
above this characteristic low-temperature scale T^ax, lattice 
and Kondo disorder scattering become indistinguishable. For 
temperatures T < T^ax, an increase in the resistivity with de- 



creasing concentration c is observed; Fig. 10(a) shows a tran- 
sition from low resistivity to an insulating behavior for T ^ 
and c ^ 0. The disorder introduced in the conduction band 
destroys the lattice coherence of the heavy quasiparticles at 
low temperatures. Even though translational invariance is re- 
stored for c ^ 0, the residual resistivity remains increasing for 
c < 0.5 due to the crossover to an insulating behavior. 

The thermoelectric power does not show much variation 
with the concentration c as can be seen in Fig. 10(b). In con- 
trast to the resistivity, the thermoelectric power depends on the 
asymmetry of the spectrum [see Eqs. (42) and (43b)]. It turns 
out that the asymmetry of the spectra as well as the asymmetry 
of the relaxation time as defined in Eq. (44) change only very 
weakly with the concentration, even though there is a signifi- 
cant renormalization of the effective media A^/g(z) [Eq. (21)] 
which enters the effective site calculations. 

Next, we investigate the transport properties as a function 
of disorder in the case in which ligand substitution keeps the 
system in a metallic regime for all values of c. In Fig. 1 1(a), 
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the resistivity is displayed for e'* 



-6ro and e? 



-4ro. 



In this case, we have a disordered conduction band 



1200 



where both subsystems A and B are asymmetric. By choosing 
a filling of «tot =1.8 — 0.2c and e4 — e? = 0.2ro, a nearly c- 
independent / occupation n/ sa 0.92 is achieved, significantly 
departed from the Kondo insulator regime. This is clearly vis- 
ible in the figure: in the limits c ~Q and c = 1, the resid- 
ual resistivity p(r^O) is small as expected for metallic sys- 
tems. For < c < 1, the disorder in the conduction band leads 
to an enhanced resistivity. The residual resistivity peaks be- 
tween c = 0.4 and 0.6. At high temperatures, the resistivity 
becomes independent of the concentration c: the incoherent 
Kondo scattering dominates the scattering processes over the 
lattice disorder. 

In Fig. 11(b), the thermoelectric power is shown for the 
same parameter set. The overall appearance is equal to the 
case with Kondo holes, but, in addition, we observe a sign 
change at small temperatures. This sign change originates 
from the change of asymmetry in the relaxation time t{(o) 
as defined in Eq. (44). To illustrate this point, we plot the an- 
tisymmetrized relaxation time F{(i)) := T(a)) — t(— o) in the 
vicinity of the chemical potential for T\ = 0.48 K as inset in 
Fig. 1 1(b). At this temperature, S{T,c = 0.1) reaches its low- 
temperature minimum, while for c = 1, a positive thermoelec- 
tric power is found. For intermediate c = 0.4, F{co) almost 
vanishes as can be seen in the inset; simultaneously, the value 
of the thermoelectric power at Ti is very small. For concen- 
trations c < 0.4, the positive slope of F{(o) leads to a negative 
thermoelectric power, while for c > 0.4, a negative slope is 
found corresponding to a positive thermoelectric power. This 
change of asymmetry is a consequence of the very subtle re- 
distribution of spectral weight also seen in the single particle 
spectra (not depicted here). 




FIG. 12: (Color online) (a) Resistivity p{T) and (b) thermoelectric 
power S(r) as functions of T for different concentrations c of system 
A; all parameters as in Fig. II, except for V^ = BFq and S/Fq = 
10-4. 



2. Resistivity and thermopower for Va ^ Vb 

In a real material, it is more likely that the hybridization 
changes with the local environment. Two different scenarios 
could be discussed. The first one considers each individual 
local configuration to be characterized by a probability distri- 
bution function for the hybridization P{V), which will lead to 
a distribution of Kondo scales. This scenario was suggested 
as a possible route to non-Fermi-liquid behavior in HFSs."'^'* 
We, however, restrict ourselves to A/B ligand disorder, yield- 
ing one effective hybridization for each type of lattice sites A 
and B: Va and Vb- For a homogeneous medium F as defined in 
Eq. (20), this leads to two different energy scales Ta and Tb for 
the local effective sites A and B. We expect a crossover from 
one heavy fermion with a lower characteristic energy scale to 
another with a higher value of Tq. 

We take Fig. 11 as reference and adapt the choice of pa- 
rameters to it. There are two possibilities to achieve Va 7^ Vb'- 
we can change Va or Vb with respect to Fig. 11. Both cases 
are discussed in the following. We reduced the small imag- 
inary part S/Fq to 10^4 in order to reduce the error in the 
residual resistivity. In Fig. 12, the resistivity and the thermo- 
electric power for ligand disorder with changing hybridiza- 



tion from a low value of V^ 
ofVi- 



3rg at c 



1 to a high value 
Vq ~ 5.64rQ are shown; in Fig. 13, the hybridization 
conversely changes from a high value of V^ = Vq ~ 5.64rQ at 
c = 1 to a low value of Vg = SFq. All other parameters are the 
same as in Fig. 11. 

In both figures, the CPA end points c = and c = 1 nicely 
illustrate the different energy scales, which are proportional 
to the temperature value of the resistivity maximum and the 
corresponding maximum of the thermoelectric power With 
varying concentration, the thermoelectric power reflects the 
crossover from one energy scale to the other by a decreasing 
height of one maximum and an increasing height of the other 
one. At intermediate concentrations, two maxima in ^(r) are 
found, which clearly shows the presence of two energy scales. 
This feature is most pronounced in Fig. 12 for c = 0.4, where 
a shallow double maximum structure is observed. 

The behavior of the resistivity is completely different in the 
two cases of Figs. 12 and 13. In Fig. 12, the residual resis- 
tivity increases with decreasing c and peaks around c = 0.3. 
Due to the change of particle numbers with the concentration, 
the residual resistivity is not symmetric around c — 0.5. At 
intermediate concentrations, we find an extended crossover 
regime. At c = 0.2, a non-Fermi-liquid behavior of the resis- 
tivity is observed over more than three decades in temperature. 
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FIG. 13: (Color online) (a) Resistivity p{T) and (b) thermoelectric 
power S{T) as functions of T for different concentrations c of system 
A; all parameters as in Fig. 11, except for V^ = SFq and S/Fg = 
10-4. 



while the thermoelectric power is already dominated from the 
B site transport. Although the c = curves in Fig. 13 have the 
same hybridization as the c = 1 curves in Fig. 12, they differ 
in the number of electrons per unit cell, since in both cases «tot 
is given by Wtot = 1 -8 — 0.2c. The change in the resistivity as a 
function of concentration appears to be much more gradually 
with a continuous and monotonic shift of the position of the 
resistivity maximum. 



VI. CONCLUSION AND OUTLOOK 

We have presented a detailed analysis of the transport prop- 
erties of disordered heavy-fermion compounds based on a 
combination of CPA to treat the disorder and DMFT to handle 
the correlation effects. Our approach reduces to the standard 
DMFT at those points corresponding to concentrated systems 
(c = and c = 1); i.e., our treatment of the transport proper- 
ties interpolates between two concentrated systems consisting 
only of pure A or B sites. 

To make contact to typical experimental situations, we 
have studied two types of local disorder: introduction of 
Kondo holes and disorder on the ligand sites. In the case 
of doping with Kondo holes, we find an increasing charac- 



teristic low-temperature scale To, defined via the quasiparti- 
cle renormalization factor, as a function of the hole concen- 
tration for constant hybridization. This is consistent with a 
previous comparison between the DMFT and SIAM energy 
scales. '^'^4 However, in experiments, often only a weak de- 
pendence of the low-energy scale as a function of doping has 
been observed. ^'^ Therefore, we investigated the change of 
the transport properties as a function of concentration for a 
fixed effective low-energy scale. Such an ansatz is consistent 
with an experimentally found increase of the lattice constants 
upon Ce substitution by La, justifying a linear reduction of 
the hybridization strength with increasing concentration of the 
Kondo holes. 

The calculated resistivity normalized to the Ce concentra- 
tion shows a characteristic logarithmic upturn when lowering 
the temperature. In the periodic system, the resistivity de- 
creases again after a maximum is reached, which is a sign of 
lattice coherence. Introducing Kondo holes reduces the scat- 
tering of the / electrons and gradually destroys lattice coher- 
ence. The normalized resistivity continuously varies from that 
of a concentrated system to that of a typical dilute Kondo scat- 
terer While we observe scaling with the concentration at high 
temperature, multiscattering of conduction-band electrons in- 
troduces an additional enhancement factor 1/(1 — ac) in the 
residual resistivity for c ^ 0, which can become already im- 
portant for concentrations less than 10%. Since the thermo- 
electric power is given by a ratio of two transport integrals, it 
is much less sensitive to disorder. 

The scaling properties of doped and undoped heavy- 
fermion materials appear to be inconsistent with a simple 
Fermi-liquid theory and have caused a lot of attention in re- 
cent years. The terminology "non-Fermi liquid" was coined 
for those materials whose thermodynamic and transport prop- 
erties are not well understood — for a review, see Ref . 1 . Usu- 
ally, the deviation from a T^ law has been attributed to scatter- 
ing of conduction electrons on thermal and quantum fluctua- 
tions of the magnetic order parameter in the vicinity of a quan- 
tum critical point. Such effects would generate a A:-dependent 
self-energy which is not included in our treatment. The ef- 
fective site for the Kondo-hole disorder problem reaches the 
strong-coupling fixed point at low enough temperature in our 
CPA-NRG theory. Consequently, we know analytically that 
the local self-energy lA (z) has Fermi-liquid properties in this 
regime. The scaling regime for such a T^ law, however, turns 
out to be very small. The self-consistency condition and mul- 
tiscattering contributions to transport properties reorganize 
spectral weight in our many-body calculation. Therefore, the 
crossover regime starts relatively early. In this regime, one 
might be able to identify a rather extended temperature range 
in which a power law with a fixed exponent a <2 can be ob- 
tained with the same accuracy as determined in experiments. 
However, the physical meaning of such an exponent obtained 
in a crossover from a low-temperature scaling regime is not 
apparent to us. Thus, we did not intend to extract such a pa- 
rameter 

In the case of ligand disorder, we have shown calculations 
which interpolate between two metallic systems as well as 
describe the crossover from a Kondo insulator to a metallic 
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HFS. Each unit cell contains one Kondo scatterer which ex- 
periences a different environment. A changing hybridization 
matrix element leads to two different local low-temperature 
scales. For a constant hybridization matrix element, we do not 
observe a shift in the position of the resistivity maximum. The 
resistivity follows the typical heavy-fermion behavior with 
a concentration-dependent residual resisitivity which peaks 
close to the concentration of c = 0.5. In the thermoelectric 
power, however, we note a sign change at very low tempera- 
tures from positive values for c < 0.5 to negative values for 
c = 0.9 similar to the one observed'^ in CeCu2Si2. The ther- 
moelectric power is very sensitive to particle-hole asymme- 
tries in the spectral functions. 

Our results clearly demonstrate that the combination of 
CPA and DMFT is capable of accurately capturing the inter- 
play of effects that are connected to local correlations, on the 
one hand, and the average influence of disorder, on the other 
hand. We discuss ligand disorder in addition to the conven- 
tionally investigated introduction of Kondo holes. Such lig- 
and disorder is relevant for a variety of HFSs and, in particu- 
lar, interesting in connection with CeCug-;, Au^, the paradigm 
of non-Fermi-liquid behavior in the vicinity of a quantum 
phase transition. Due to the lack of nonlocal fluctuations in 



the DMFT-CPA treatment, it will not be possible to study the 
quantum critical regime within our method. However, within 
our approachO we can investigate the influence of the lig- 
and disorder on magnetic properties to obtain the magnetic 
phase diagram. Another aspect is the role of orbital degrees of 
freedom, in particular, on thermoelectric properties of HFSs. 
Again, the method introduced here is valid for such an ex- 
tended model, too, but the calculations become much more 
involved. Investigations along these lines are in progress. 
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